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Pellet  gains  in  excess  of  100  will  probably  be  necessary  for  most  applications  of  inertial  fusion.^ 
In  order  to  achieve  these  high  gains  a  number  of  critical  physics  elements  must  be  controlled.  These 
include  (1)  high  coupling  efflciency,  (2)  low  fuel  preheat,  (3)  implosion  symmetry,  (4)  implosion 
stability  (ablation  pressure)  and  (5)  an  ignition  scheme.  These  factors  interact  with  each  other  pro¬ 
viding  conflicting  requirements.  The  purpose  of  this  paper  is  to  report  on  the  theoretical  and 
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ANALYSIS  OF  STABILITY  AND  SYMMETRY  IMPLICATIONS  FOR  ICF 


I.  INTRODUCTION 

Bsllet  gains  in  excess  of  100  will  probably  be  necessary  for 
most  applications  of  inertial  fusion.*  In  order  to  achieve  these 
high  gains  a  number  of  critical  physics  elements  mist  be  control¬ 
led.  These  include  (l)  high  coupling  efficiency,  (2)  low  fuel 
preheat,  (3)  inplosion  symmetry,  (U)  inplosion  stability  (ablation 
pressure)  and  (5)  an  ignition  scheme.  These  factors  interact  with 
each  other  providing  conflicting  requirements.  In  particular  the 
first  two  items  are  directly  in  conflict  with  the  second  two.  For 
exanple,  high  coupling  efficiency  and  low  fuel  preheat  requires 
control  of  deleterious  plasma  instabilities.  These  instabilities 
generally  scale  as  IX^.  Therefore  they  are  usually  controlled  by 
the  use  of  lower  intensities  or  shorter  laser  wavelength.  But 
lower  laser  intensities  are  generally  associated  with  thinner  shell 
or  double  shell  target  designs,^  and  these  higher  aspect  ratio 
designs  place  severe  requirements  on  laser  symmetry  and  target 
stability.  Smoothing  out  laser  nonuniformities  by  lateral  thermal 
conduction^  requires  that  the  absorption-to-ablation  distance  be  on 
the  order  of  the  target  radius .leading  to  longer,  not  shorter  wave¬ 
lengths.  This  separation  distance  also  produces  a  lower  hydro- 
dynamic  efficiency.  The  control  of  the  Rayleigh-Taylor  instability 
by  use  of  smaller  aspect  ratio  targets  increases  the  ablation  pres¬ 
sure  requirements,  and  therefore  the  laser  intensity,  thus  reintro¬ 
ducing  the  preheat  problem.  iJse  of  shorter  wavelength  lasers  to 
reduce  the  plasma  instabilities  also  reduces  the  separation 
distance  between  the  absorption  and  ablation  surfaces  thus  reducing 
the  effect  of  thermal  smoothing.  The  problem  is  then  to  find  a 
cortpromise  among  the  conpeting  effects  which  will  lead  to  an 
acceptable  high  gain  target.  The  puirpose  of  this  paper  is  to 
report  on  the  theoretical  and  conputational  work  that  has  been 
carried  out  in  connection  with  the  experimental  work  at  NRL  in 
evaluating  the  design  constraints  that  are  inposed  primarily  by  the 
symmetry  and  stability  requirements  on  directly -driven  laser  fusion 
targets. 
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We  have  approached  this  problem  with  the  viewpoint  that  suf¬ 
ficiently  low  intensities  will  be  required  that  classical  physical 
processes  will  dominate.  Thus  we  avoid  as  nuch  as  possible  the 
problems  of  fast  electron  production  and  nonclassical  ener^  trans¬ 
port  mechanisms.  Our  philosophy  is  then  to  find  what  limits  the 
critical  elements  will  place  on  the  physically  relevant  parameters 
such  as  target  radius,  aspect  ratio,  laser  intensity  and  laser 
wavelength.  We  have  developed  scaling  laws  that  suggest  that  there 
is  a  possible  window  near  10^“*  W/cm^  and  near  lum  laser  wave¬ 
lengths.**  In  deriving  our  scaling  laws  we  have  implicitly  assumed 
that  what  we  desire  is  a  maximum  pellet  gain  Q,  and  that  any  laser 
energr  necessary  to  drive  the  target  can  be  achieved  by  building  a 
sufficiently  large  laser.  laser  energy  scales  as  R^(dR/R)  and  is 
strongly  dependent  on  target  size  but  only  weakly  dependent  on 
irradiance.  The  constraints  that  we  consider  are  limits  itqposed  by 
symmetry,  stability,  coupling  efficiency,  preheat,  and  implosion 
velocity.  We  have  not  yet  investigated  the  constraint  inposed  by 
an  ignition  design  as  an  acceptible  design  yet  eludes  us. 

Our  results  are  primarily  based  on  numerical  sinulations  with 
many  of  the  results  verified  by  experiments  performed  at  NRL. 
Detailed  conparisons  of  code  results  with  experiments  are  presented 
elsewhere  in  this  proceedings.^  We  enploy  several  different  types 
of  codes.  For  the  symmetry  and  stability  issues,  we  use  an  adap¬ 
tively  gridded,  two-dimensional,  fully  nonlinear,  Eulerian  computer 
model  (FAST2D)  and  a  triangular  grid  Lagrangian  conputer  model 
(SPLISH).  Laser  matter  coupling  and  scaling  laws,  relating  the 
asymmetry  and  stability  results  to  spherical  targets,  are  studied 
using  a  one-dimensional,  sliding  zone,  adaptively  gridded,  fluid 
code  (FASTID).  Principle  results  to  be  reported  here  include: 

(1)  requirements  for  thermal  smoothing  of  laser  nonuniformities. 

(2)  reduced  Payleigh-'Ih.ylor  growth  rates  and  saturation  via  the 
nonlinear  Kelvin-Helmholtz  shear  flow,  (3)  a  new  theory  based  on 
vortex  shedding  that  explains  the  reduced  linear  Rayleigh-Thylor 
growth  rates,  and  {k)  a  possible  wavelength-intensity  window  for 
high  gain  direct  illumination  where  the  requirements  of  target 
velocity,  symmetry,  and  now  also  efficiency  can  be  simultanesouly 
met. 


II.  MODELS  USED 

The  symmetry  and  Rayleigh -Taylor  stabilty  analysis  were  per¬ 
formed  using  the  FAST2D  laser-shell  simulation  code.  This  is  a 
fully  two-dimensional  nonlinear  Eulerian  fluid  code  ertplcying  a 
sliding  Cartesian  grid  with  variable  grid  spacing.  The  system 
enployes  1*0  zones  transverse  to  the  laser  beam  and  120  zones 
parallel  to  the  beam.  For  the  R-T  runs  we  enplcyed  a  fine  grid  of 
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0.25  um  for  ten  zones  on  either  side  of  the  ablation  layer, 
increasing  to  a  uniform  2  lun  spacing  for  most  of  the  rest  of  the 
grid.  The  code  has  been  run  vith  the  fine  grid  spacing  as  small  as 
0.10  um  with  no  noticeable  difference  in  the  results.  The  refined 
subzoning  moves  with  the  ablation  front  throu^out  the  course  of  a 
run.  The  zones  in  the  underdense  plasma  beyond  the  critical  sur¬ 
face  and  in  the  low-density  rear  portion  of  the  foil  are  exponen¬ 
tially  stretched  so  as  to  minimize  the  influence  of  the  boundaries. 
The  code  nay  be  run  in  either  planar  or  cylinderiral  coordinates 
with  either  periodic  or  outflow  boundary  conditions  in  the  trans¬ 
verse  direction.  An  option  is  also  a-zailable  to  stretch  the  outer 
zones  in  the  transverse  direction. 

FAST2D  solves  the  time -dependent  equations  of  motion  for  the 
continuity  of  mass,  momentum  and  total  energy  for  a  fully  ionized 
ideal  gas  using  the  time  split,  flux-corrected  transport  algo¬ 
rithms.  Thermal  conduction  is  incorporated  in  a  single  temperature 
model  using  an  irolicit  solution  of  the  two-dimensional  classical 
nonlinear  (y  thermal  conduction  equation  with  the  standard 

2raginskil  coefficient.  Initial  conditions  are  generated  from  a 
one-dimensional,  analytic,  quasi-atatic  equilibrium  model. °  The 
analytic  model  provide  a  solution  which  changes  only  about  .5*  when 
FAST2D  is  tested  in  one -dimensional  nearly -steady -state  mode. 

Interior  layer  Bayleigh-Thylor  growth  is  investigated  using  a 
Lagrangian  fluid  ctrnamics  model  (SPLISH)  which  bases  its  fluid 
representation  on  a  two-dimensional,  dynamically  reconnectable  grid 
of  triangles.  A  unique  reconnection  algorithm  allows  additional 
triangles  to  be  added  along  the  shear  layer  as  it  stretches  in  the 
nonlinear  phases  of  the  Hayleigh-Taylor  growth  and  enables  the  roll 
up  to  be  followed  far  into  the  nonlinear  phase  where  droplets  break 
off  and  mixing  occurs.  This  model  employs  ideal  nonviscous, 
nonconducting  hydrodynamics  appropriate  to  the  high  density 
interior  flow. 

We  investigate  the  laser  plasma  coupling  and  scaling  phenomena 
using  our  one-dimensional,  adaptively  gridded,  spherical  hydro- 
i^namics  code  (FASTID).  The  time -dependent  equations  for  continu¬ 
ity  of  mass,  momentum,  and  total  energy  are  integrated  numerically 
with  sliding  zone,  Eulerian  flux-corrected  transport.  A  flux- 
conser'/ing  scheme  is  used  such  that  the  fluxes  between  zones 
results  from  the  difference  in  fluid  and  grid  velocities.  The 
outer  zones  and  material  interfaces  move  in  a  Lagrangian  manner 
while  intermediary  zones  move  in  a  prescribed  smooth  continuous 
manner  so  as  to  maintain  the  desired  resolution  near  the  ablation 
and  critical  surfaces.  !fear  the  ablation  surface  the  zone  size  is 


held  fixed.  From  the  ablation  surface  to  the  critical  surface  the 
mass  per  zone  is  held  nearly  fixed  resulting  in  increasing  zone 
size.  Beyond  the  critical  surface  the  mass  per  zone  is  allowed  to 
increase  exponentially.  This  technique  allows  a  grid  resolution  of 
a  few  tenths  of  a  micron  near  the  ablation  surface  and  a  smooth 
transition  to  larger  zones  in  the  underdense  plasma.  It  has  a 
distinct  ad-/antage  over  Lagrangian  codes  which  find  it  difficult  to 
enhance  the  resolution  where  needed  without  wasting  zones  else¬ 
where. 


Thermal  conduction  is  computed  using  the  usual  reciprical 
average  of  the  classical  and  limiting  heat  flux,  1/F  =  l/F^^+l/F^^ 
where  F^^  =  f{kT/m^)  ^'^^n^kTVT/I  7T!  .This  flux  limiter  F<-^  is  a 
construct  commot.ly  used  in  fluid  codes,  for  which  there  is  not  as 
yet  experimental  verification  for  reactor  size  plasmas.  Ve  have 
used  values  of  f  ranging  from  0.6  to  0.03.  Our  results  are  mainly 
quoted  for  the  Bell-Evans-Nicholas^  value  of  f  =  0.1.  The  results 
are  insensitive  to  values  of  f  >  0.1  as  long  as  the  ;4ach  number  at 
the  critical  surface  is  also  not  too  much  greater  than  one  (small 
spherical  divergence).  laser  light  is  absorbed  by  inverse  brems- 
strahlung  of  both  incoming  and  reflected  light.  The  fraction  of 
absorbed  light  not  deposited  by  the  inverse  bremsst rah lung  is 
deposited  smoothly  over  two  zones  including  the  critical  surface  so 
as  not  to  generate  artificial  shocks. 

The  equations  are  integrated  in  strict  conservation  form. 
Pressure  and  temperature  are  thus  derived  by  subtract'ng  directed 
kinetic  energy  from  total  energr.  The  fluid  equations  are  closed 
with  an  equation  of  state  expressing  ?  =  F(F,p)  and  T  =  G(E,p). 

Mear  the  ablation  surface  these  equations  nust  include  the  effects 
of  solid  material,  Fermi -temperature,  and  fractional  ionization. 

We  base  our  equation  of  state  on  the  analytic  model  (AHEOS)  in  the 
THARTD  code,®  and  expansions  of  the  Fermi -Dirac  functions  by 
Clayton.  This  model  goes  smoothly  from  a  Thomas -Fermi -Dirac  model 
for  degenerate  electrons  at  high  density  and  low  temperature  to  a 
Saha  equilibrium  at  high  temperatures  and  low  densities.  Such  a 
treatment  is  necessary  in  order  to  correctly  calculate  the  acceler¬ 
ation  of  near-solid,  low-teiperature  target  shells  without  having 
the  shell  compress  or  explode  unrealistically. 

ITT.  SYMMETRY 


One  of  the  critical  elements  for  high  gain  target  designs  is 
the  degree  of  symmetry  that  must  be  maintained  in  the  implosion 
process.  For  high  gain,  the  asymmetry  of  the  ablation  pressure 
must  be  better  than  a  few  percent  (6P/?  <  1-2*).^  Either  the  laser 
illumination  mst  be  uniform  to  this  percentage,  or  some  additional 
smoothing  mechanism  is  required.  As  proposed  by  Nuckolls,®  lateral 
thermal  conduction  could  smooth  out  the  laser  nonuniformities 
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Figure  1.  The  •/aria':ion  in  ablation  pressure  (ilp/<p>}  for  all  six 
scale  lengths  plotted  as  a  function  of  the  average 
laser  intensity  (<I>)  and  as  a  function  of  the  average 
distance  from  the  ablation  surface  to  the  critical 
surface  (<Dg^5,>). 


yielding  oP/?  <<  dl/I  provided  the  separation  of  the  absorption 
region  and  the  ablation  surface  is  sufficiently  large  ~ 

vhere  X^  is  the  asymnetry  wavelength.  We  have  performed  a 
quantitative  assessment  of  this  smoothing  effect  using  our  two- 
dimensional  laser-shell  simulation  code  FAST2D. 

Pesults  are  presented  for  a  thin  (15  ua  thick)  plastic  (CH) 
foil  accelerated  to  quasi  steady  state  (6  nanoseconds)  by  a  laser 
beam  with  an  intensity  •-’ariation  of  ±ho*  „  ~  ~  O.Sol 

and  a  laser  wavelength  of  l.Co  pra.  In  these  simulations  the 
initial  conditions  were  provided  by  a  quasi-static  equilibrium 
model  with  the  relatively -low  average  laser  intensities  of  <I>  = 
5x10^2^  1x10^^,  and  2xlC^^  W/cm^.  The  laser  beam  profile  was 
'laussian  with  a  reduced  intensity  in  the  center  of  the  beam.  Since 
we  are  only  interested  in  laser  asymmetry  effects,  we  invoke  the 
periodicity  of  the  code  and  only  model  peak  to  peak  variations. 

The  code  was  run  with  the  wavelength  of  the  laser  inhoraogeneity 
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Figure  2.  Coirparison  of  the  numerical  results  for  the  ICC- 

and  200un  scale  lengths  with  the  luC-um  scale  length 
experimental  results  (Fef.  12)  and  with  the  llO-um 
theoretical  results  'dashed  curve,  Fef.  13). 


varying  from  X,  =  100  microns  to  X^  =  600  microns.  The  results 
are  summarized  in  Fig.  1,  where  the  ’ra,riation  in  ablation  pressure 
is  plotted  as  a  ^'unction  of  the  average  laser  intensity,  or  equiva- 
Isntly  the  average  distance  between  the  ablation  and  critical  sur- 
f  faces.  The  direct  correlation  between  and  the  laser 

■  intensity  has  been  noted  both  experimentally’’^  and  theoretically^®. 


We  find  that  lateral  energy  flow  in  the  form  of  transverse 
thermal  conduction  does  smooth  the  effects  of  laser  asymmetries.^^ 
The  degree  of  smoothing  is  directly  related  to  the  distance  between 
the  ablation  and  critical  surfaces.  The  results  are  in  agreement 
with  experimental  results  obtained  at  the  ?fa.val  Research 
Laboratory^ ^  and  with  an  analytic  model  developed  at  WRL^®  as  shown 
in  Fig.  2  where  we  compare  the  foil  velocity  nonuniformity 
(V.,^^/V.^,j-1)  for  the  100  micron  and  200  micron  asymmetries  to 
experiniental  results  for  a  1^0  micron  asymmetry  with  "  2. 

In  Fig.  3a -b  we  present  a  detailed  look  at  the  case  where  <1>  = 
1x10^3  v/cm^  for  X^  *  100  microns  and  600  microns.  Pressure 
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<1  >  •  I  0  »  10  v»/cm‘ 

X^>  6  0  0  >i.  m 


"ijure  3. 


a)  Pressure  contours  for  <I>  =  V?/cn:^  and  Xg  * 

6Q0  ym  increnents  of  0.5  of  aaxinun  pressure.  laser 
comes  in  from  the  right.  Mote  also  critical  and  abla¬ 
tion  surface,  b)  Pressure  contours  for  -  100  um. 


contours  are  in  increments  of  0.1  counting  from  the  outside  inward, 
and  the  laser  light  is  coming  in  from  the  right.  For  the  600 
micron  case  severe  distortion  of  the  foil  has  taken  place.  In  the 
100  micron  case  the  ablation  surface  is  essentially  equi'/alent  to  a 
contour  surface  and  sufficient  smoothing  has  occurred  to  reduce  the 
ablation  pressure  variation  to  0.5^.  This  results  is  basically 
intuitive.  If  the  scalelengths  of  the  transverse  tenperature  vari¬ 
ations  are  comparable  to  the  separation  distance,  then  a  thermal 
wave  propagating  transverse  to  the  laser  beam  should  smooth  out 
moderate  temperature  variations  in  about  the  time  that  the  laser 
energy  is  conducted  to  the  ablation  surface.  These  simlation 
results  allow  us  to  obtain  quantitative  estimates  of  the  require¬ 
ments  for  thermal  smoothing  which  can  be  used  in  scaling 
arguments. 
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IV.  STABILITY 

The  stability  of  the  implosion  is  another  critical  element  for 
high  gain  targets.  The  use  of  lower  laser  intensities,  and  there¬ 
fore  lower  ablation  pressures,  to  reduce  the  plasma  instability 
effects  requires  the  use  of  higher  aspect  ratio  (thinner  shelled) 
targets.  We  have  investigated  the  impact  that  the  Rayleigh-Taylor 
instability  has  in  limiting  this  aspect  ratio. 

This  instability  may  grow  in  the  ablation  layer  or  at  internal 
interfaces  where  materials  of  different  densities  abut,  whenever 
the  pressure  and  density  gradients  have  opposite  sign.  Figure  7 
shows  four  snapshots  of  a  Rayleigh-Taylor  unstable  interface  calcu¬ 
lated  using  our  SPLISH  model.  An  initially  sinusoidal  perturbation 
becomes  nonlinearly  distorted  into  a  classical  "bubble  and  spike" 
pattern  as  the  exponential  growth  phase  begins  to  saturate.  When 
sufficient  shear  is  generated  in  the  flow  by  the  RT  mode,  the 
Xelvin-Helmholtz  instability  (KH)  becomes  nonlinearly  active  and 
broadens  the  spike  tip,  eventually  spreading  the  falling  spike 
material  out  and  rolling  it  up  around  the  large  vorticities  which 
were  generated  in  the  density  gradient  and  subsequently  advected 
down  the  sides  of  the  spikes. 

We  have  modeled  the  ablation  layer  Rayleigh-Taylor  problem 
with  our  two-dimensional  laser-shell  simulation  code  FAST2D. 

Initial  conditions  are  again  obtained  from  the  quasi-static  equi¬ 
librium  model.  We  have  modeled  the  RT  instability  far  into  the 
nonlinear  regime  beyond  the  point  of  foil  fragmentation  for  two 
types  of  initial  perturbations.  One  is  a  multinode  perturbation 
where  the  first  ten  modes  are  excited  with  equal  amplitudes  and 
random  phases  corresponding  to  an  initial  density  perturbation  of 

“  0.08).^**  The  other  is  single  mode  sinus¬ 
oidal  perturbations  along  the  density  peak  corresponding  to  a 
total  initial  density  perturbation  of  t  ^  ~  0.023]. 

The  wavelength  (A)  is  varied  in  the  range  ^2  V  AR  <  10  where  AR 
is  the  cold  foil  thickness  (20  microns).  In  both  multi-mode  and 
single  mode  cases  we  find  linear  growth  rates  below  classical 
values  by  a  factor  of  two  to  three,  and  a  cutoff  in  the  growth 
rates  for  wavelengths  less  than  the  cold  foil  thickness  (see 
Fig.  4). 

A  striking  result  we  found  is  the  dominance  of  nonlinear 
effects  in  determining  the  largest  viable  aspect  ratio  for  thin 
shell  targets.  Although  the  linear  growth  rates  increase  as  k^^^ 
almost  up  to  the  cutoff,  the  Kelvin-Helmholtz  (KH)  rollup  dominates 
at  large  k,  reducing  the  penetration  rate  of  the  dense  spike  below 
the  free  fall  value  and  effectively  doubling  the  aspect  ratio  of 
the  foil. 
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Figlire  A  coinparison  of  the  con^Jutational  growth  rates 

to  the  classical  value  {(kg'^''^, - )  and  to  the  vortex 

shedding  theory  (0).  Cold  foil  thickness  is  2C  urn. 


Ihe  growth  rate  of  the  perturbation  shifts  fron  exponential  to 
quadratic  (free-fall)  when  the  amplitude  of  the  disturbance  is  on 
the  order  of  1/5  the  wavelength  (a  ~  0.2  Xl .  Ihis  results  in  a 
different  saturation  level  for  the  short  and  long  wavelength  modes. 
Ihe  amplitude  of  the  large  wavelength  modes  (  X/AP.  >  1.25)  goes  into 
"free-fall"  at  a  later  time  and  with  an  effective  "free-fall" 
acceleration  comparable  to  the  foil  acceleration  fl.5>'10^^  cm/s^). 
The  exponential  growth  of  the  short  wavelength  modes  ( X/AR 
<  1.25)  saturates  earlier  in  time  (since  the  linear  growth 
rates  are  coiqparable  and  the  wavelengths  shorter)  and  the  "free- 
fall"  acceleration  of  the  spike  is  reduced  by  more  than  a  factor  of 
two  over  the  long  wavelength  nodes.  Ihis  is  evident  from  Fig.  5 
where  both  the  "free-fall"  acceleration  of  the  spike  and  the  final 
aspect  ratio  of  the  foil  are  plotted  as  a  function  of  perturbations 
wavelength  (X/AR)  for  single  mode  perturbations.  This  reduced  rate 
of  penetration  of  the  dense  spike  below  its  "free-fall"  value  is  a 
direct  consequence  of  the  tip-widening  brought  on  by  the  ?CH  rollup. 
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Fl^re  5 


A  plot  of  the  "free-fall"  acceleration  ("g"  in  units  o 
of  the  spike  and  the  final  aspect  ratio  of 


cn/s  2 ' 


the  foil  as  a  function  of  the  perturbation  vavelength 


(X/AR).  The  solid  circles  are  the  aspect  ratios 
obtained  with  a  laser  intensity  of  5*10^^  W/cm^  and  th 
^  synbol  is  the  aspect  ratio  obtained  with  a  perturba¬ 
tion  wavelength  of  200  ym. 


This  nonlinear  result  alters  the  conclusions  one  draws  with 
respect  to  shell  integrity  as  a  function  of  perturbation  wave¬ 
length.  As  a  result  of  the  earlier  saturation  of  the  exponential 
growth,  and  the  reduced  "f ree-fall"growth,  the  aspect  ratio  (R/AR) 
allowable  to  the  short  wavelength  perturbations  has  been  increased 
by  a  factor  of  two  over  the  long  wavelength  todes.  AR  is  the  cold 
foil  thickness  (20  ura)  and  R  is  the  distance  traveled  until  the 
foil  uniformity  is  destroyed.  This  is  defined  as  the  point  at 
which  either  35%  of  the  foil  mss  has  been  ablated  or  the  rear  of 
the  foil  was  no  longer  accelerated  uniformly.  The  solid  circles 
are  the  aspect  ratio  obtained  with  a  laser  intensity  of  5*0x10^^ 
W/cn^,  where  the  same  effect  is  evident.  The  #  symbol  is  the 
aspect  ratio  obtained  for  a  200  um  perturbation  (1=1.0*10^^  W/cm^). 
As  one  would  intuitively  guess,  the  allowable  aspect  ratio 
increases  again  for  perturbation  wavelengths  very  nuch  larger  than 
the  foil  thickness.  Although  the  linear  growth  rate  increases  as 
k^^2  up  to  the  cutoff,  nonlinear  effects  in  the  form  of  the  KH 
roll-up  are  quite  effective  in  saturating  the  RT  instability. 


reducing  the  rate  of  "fall"  of  the  spike,  and  increasing  the  life¬ 
time  of  the  foil. 

The  vorticity,  which  plays  such  a  dominate  role  in  the  non¬ 
linear  evolution  of  the  RT  instability,  is  also  responsibile  for 
the  reduced  linear  growth  rates.  In  particular,  the  shedding  of 
the  vorticity  generated  at  the  ablation  layer  is  similar  to  the 
shedding  of  vorticies  that  arise  from  the  flow  about  a  bluff  body. 
Not  only  does  this  shedding  reduce  the  linear  RT  growth  rate  but  a 
sj-mnetric  von  Karman  vortex  street  develops  in  the  blow  off  region. 
The  symmetric  vortex  street  is  very  unstable  and  thus  could  also 
have  serious  effects  on  laser  absorption,  thermal  conduction  and 
instabilities  in  the  underdense  region. 

Figure  6  shows  the  vortex  shedding  process  from  the  ablation 
layer  for  the  20  um  perturbation  wavelength  when  the  system  is  in 
the  linear  regime.  The  finite-sized  vortices  are  generated  by  the 
baroclinic  (non-colinear  density  and  pressure  gradients)  nature  of 
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Figure  6.  Vorticity  contours  for  the  20  um  mode  (X/AR  =  1)  when 
the  system  is  in  the  linear  regime.  The  contours  are 
in  10%  increments  of  the  maximm  vorticity  (+2.75>‘10®, 

-  2.87x10®  S“M.  The  solid  (dashed)  contours  are  posi¬ 
tive  (negative)  vorticity.  The  heavy  lines  delineate 
the  10%  (right)  and  90%  (left)  density  contour. 


the  profiles.  Uie  vorticity  generated  at  the  ablation  surface  is 
then  convected  with  the  flow  and  shed  into  the  blowoff.  This  shed¬ 
ding  process  has  a  measured  Strouhal  number  (dimensionless  shedding 
frequency)  which  agrees  very  well  with  theory.  As  a  result  of  this 
shedding,  the  classical  RT  growth  rate  is  reduced.  We  account  for 
this  reduced  growth  rate  by  writing  y*  =  (kg)^^2-(xf,  where  f  is  the 
measured  shedding  frequency  and  a  is  some  fraction  (less  than  one). 
We  choose  a  *  0,56  to  fit  the  20  un  mode.  The  coii5>utational 
results  (— •— )  are  cot!5)ared  with  the  classical  growth  rate 

((kg)^''2^  - )  new  theoretical  results  (o).  The  agreement 

between  the  computatioral  results  and  our  vortex  shedding  theory'-  is 
excellent  (see  Fig.  U). 

Vorticity  continues  to  be  shed  from  the  ablation  layer  until 
the  spike  amplitude  becomes  con^jarable  to  half  the  perturbation 
wavelength.  At  this  point  the  bubble  region  ceases  to  ablate  and 
the  vorticity  is  advected  down  the  sides  of  the  spikes  where  it 
sets  up  a  reversed  flew.  Tt  is  this  secondary  flow,  the  '<H  rollup, 
which  causes  the  tips  of  the  spikes  to  widen  and  reduces  their  fall 
rate.  It  is  apparent  that  vortex  dynamics  is  playing  a  dominant 
role  in  both  the  linear  and  nonlinear  phases  of  the  RT 
instability. 

The  interior-layer  Rayleigh -Taylor  problem  can  occur  at  many 
different  times  and  locations  during  the  implosion  of  a  complicated 
nultilayer  target.  The  preignition  mix  problem  is  the  most 
familiar  situation.  This  occurs  near  the  final  stages  of  compres¬ 
sion  of  the  fuel  when  the  interface  between  the  higher  density 
pusher  and  the  fuel  compress ionallj'  decelerates.  A  three-layer  mix 
problem  arises  where  a  light  fuel  layer  is  separated  from  a  denser 
pusher  shell  by  a  thin  but  very  dense  heat-shield  layer.  Short- 
wavelength  perturbations  permit  the  upper  and  lower  surfaces  of  the 
heat  shield  to  move  independently,  so  the  RT  instability  can  result 
even  though  overall  stability  is  expected  at  long  wavelengths.  In 
order  to  investigate  the  nonlinear  regime  of  the  RT  instability  for 
these  two-  and  three-layer  mix  problems  we  have  developed  a 
lAgrangian  model  (SPLISH)  which  bases  its  fluid  representation  on  a 
two-dimensional  dynamically  reconnectable  grid  of  triangles^ 

In  the  three  layer  mix  problem  the  densities  of  the  layers  are 
in  the  ratio  G. 1:2:1  so  that  the  unstable  interface  has  an  Atwood 
number  of  1/3.  At  this  Atwood  number  the  evolution  of  the  unstable 
interface  is  governed  by  the  interaction  between  the  RT  and  KH 
instabilties  in  both  the  linear  and  nonlinear  regimes.  The  evolu¬ 
tion  of  the  RT  and  KH  instabilities  can  be  separated  into  four 
distinct  phases:  (l)  the  initial  classical,  linear  growth  phase, 
ending  when  the  spike-to-bubble  height  is  half  the  perturbation 
wavelength,  or  for  thin  layers,  when  the  spike  to  bubble  height 
equals  the  layer  thickness;  (2)  a  saturated  RT  growth  phase  with 
steadily  increasing  KH  rollup  and  increasing  wavelength,  ending 
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Figure  T.  Four  stages  of  the  nonlinear  evolution  of  a  three  layer 
mix  P.ayleigh-Taylor  instability  problem.  Density 
ratios  are  0.1:?:1. 

when  the  spike  tip  becomes  wider  than  the  spike  body;  (3)  a  small 
scale  structure  phase;  (It)  a  turbulent  mixing  phase  (see  Fig.  7). 
The  third  phase  is  characterized  by  downward  flow  from  the  thinned 
layers  above  the  bubbles  resembling  Jets  feeding  the  spike.  A 
stagnation  point  forms  at  the  point  where  the  Jets  collide,  forcing 
an  upward  Jet  opposed  to  the  spike  fall.  In  the  last  phase,  large 
amplitude  waves  run  across  the  stable  interface  as  well,  giving 
rise  to  Jetting  and  wave-breaking.  Many  similarities  exist  between 
the  ablative  and  nonablative  cases.  The  KH  instability  plays  a 
major  role  in  both  cases  in  that  it  reduces  the  penetration  rate  of 
the  dense  spike  below  its  free-fall  value. 

V.  SCALINO  UWS 


In  order  to  relate  the  symmetry  and  stability  results  to  the 
design  of  large  reactor  size  targets,  we  require  scaling  laws  rela¬ 
ting  the  laser  parameters  to  the  target  parameters.  We  have  used 
our  one-dimensional  spherical  hydrof^namics  code  to  determine  the 
scaling  laws  for  (a)  separation  distance  of  ablation  to  critical 
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surface,  and  (b)  ablation  pressure  as  a  function  of  laser  intensity 
and  wavelength.**  Rather  than  attenpt  to  design  a  set  of  optimized 
implosions  in  these  scaling  runs  we  used  pure  CH  shells  of  constant 
thickness  125  uni»  We  used  no  prepulses  to  generate  artificially 
large  initial  separation  distances  and  we  ran  each  case  until  a 
quasi-steac^  state  was  established.  The  results  are  expressed  in 
terms  of  absorbed  irradiance  since  the  absorption  process  is  not 
well  understood  for  reactor  scale  plasmas.  In  our  simulations  we 
used  pellet  radii  of  1,2,  and  5  millimeters.  Since  laser  energy 
scales  as  R^{AR/R)  a  much  larger  radii  would  imply  either  excessive 
laser  ener^  or  aspect  ratio.  Smaller  targets  would  not  contain 
sufficient  fuel  for  efficient  burn  and  high  gain.  'Ven  a  1  mm 
radius  pellet  is  marginally  small.  laser  wavelengths  of  0.53,  1.06 
and  2.7  microns  in  the  range  of  current  laser  technology  were 
chosen,  and  intensities  were  varied  from  10^ ^-10^^  V/cm^.  Figure  8 
shows  a  least  squares  fit  to  the  ablation  pressure,  as  a  function 
of  laser  intensity  and  wavelengt.h  and  target  radius  of  the  form 
P~(lg^X'J®«25R-0,18),7  Figure  9  shows  D/R  as  a  function  of  laser 
intensity  for  three  wavelengths  and  two  radii.  This  data  can  be 
represented  by  a  least  squares  fit  of  the  form  D/R 

iif. 


(I„  aJ*8)0*7/(R+B).  In  the  region  of  interest  B  ■  1 


mliLimeter  and  is  only  weakly'  dependent  on  I^  and  Aj^. 
form  has  the  advantage  of  going  to  the  correct  limits  for  both 
large  and  small  radii.  This  surprisingly  strong  '’avelength  scaling 


This 


ABSORBED  IRRADIANCE  l.t,.  (W/CM^I 

Figure  o.  Ablation  surface  pressure  scaling  for  several  laser 
wavelengths  and  target  radii. 
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A  =  1.06 
A  =  0.53 


A  =  0.53 
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Figure  ?.  Patio  of  distance  from  ablation  to  critical  surfaces 

for  several  laser  wavelengths  and  target  radii.  A  flux 
limit  of  f=0.03  reduces  D/H  by  roughly  a  factor  of 

is  the  result  of  two  factors:  longer  wavelength  light  is  absorbed 
at  lower  densities,  and  these  lower  densities  have  longer  scale 
lengths.  Sven  stronger  wavelength  scaling  has  been  found  in  the 
case  of  olanar  ablation.^ ^ 


Use  of  a  strong  flux  inhibitor  only  weakly  affects  the  abla¬ 
tion  pressure  but  strongly  effects  the  separation  distance  from 
ablation  to  critical  surface.  For  instance,  use  of  a  flux  limit 
f=0.03  vcommonly  quoted)  decreases  this  separation  distance  by  a 
factor  of  3-4. 

The  hydro dyna-mic  inplosion  efficiency  is  another  critical  ele¬ 
ment  for  high  gain  systems.  In  performing  these  simulations  we 
varied  the  shell  ti^ickness  and  held  the  final  inplosion  velocity 
fixed  at  V^=2.2xl0  cra/sec.  If  most  of  the  kinetic  energy  can 
be  efficiently  converted  to  conpressional  energy,  then  this  would 
be  sufficient  ener®r  to  compress  the  taget  to  1000  times  solid 
density  at  twice  the  Fermi -degenerate  ener®^.  If  a  relatively 
unshaped  pulse  of  fixed  wave  length  is  used,  then  the  hydrodynamic 
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LASER  WAVELENGTH  (mn) 

?igure  n.  Wave length -intensity  boundaries  determined  by  the 

symmetry  requirements  (downward  sloping  solid  curves) 
and  the  stability  limits  (aspect  ratio  <  20)  for  two 
radii  pellets  (.1  and  .5  cm).  Constant  efficiency  of 
5  and  are  shown  by  dashed  line.  Operating  region  is 
above  and  to  the  right  of  the  stability  and  symmetry 
boundaries. 

efficiency  becomes  very  poor  as  the  target  radius  approaches  col¬ 
lapse.  For  this  reason  we  choose  our  aspect  ratio  such  that  the 
find  velocity  is  achieved  after  the  shell  has  inploded  to  of  its 
initial  radius.  The  laser  driver  is  then  turned  off.  We  then 
varied  the  aspect  ratio  with  intensity  and  wavelength  to  find  the 
hydrodynamic  icplosion  efficiency  given  these  constraints.  The 
efficiency  was  found  to  be  relatively  insensitive  to  the  target 
radius  for  a  fixed  aspect  ratio.  However,  for  fixe.>  separation 
distance  from  ablation  to  critical  surface,  the  efficiency 
decreased  for  decreasing  laser  wavelength!  This  result  is  the 
opposite  of  what  some  of  our  colleagues  have  predicted.  It 
results  because  shorter  wavelengths  require  higher  laser  intensity 
(given  a  fixed  smoothing  requirement).  Higher  laser  intensities 
produce  higher  blowoff  velocities,  and  thus  lower  hydrodynamic 
efficiency.  Figure  10  sumarizes  the  results  of  the  efficiency 
study  ty  showing  contours  of  equal  efficiency  superinposed  on  plots 
of  equal  separation  ratio  and  constant  ablation  pressure. 
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VI.  IMPLICATIONS  FOR  ICF 


High  gain  targets  vill  require  inplosion  velocity  unifornity 
to  better  than  1-2  percent.^  Our  flat  geometry  symmetry  results 
show  thao  if  the  separation  distance  D  is  more  then  0.3  times  the 
laser  inhonogeneity  wavelength,  then  5?/P<0.1 6I/I,  providing  suf¬ 
ficient  smoothing  for  ±10*  laser  nonuniformities  (Fig.  l).  We 
expect  roughly  the  same  factor  of  ten  smoothing  in  spherical 
geometry  if  D/R>0.3.  Ihus  one  boundary  of  our  design  region  is 
given  by  D/?.sC.3  as  shown  in  Fig.  10  for  the  two  limits  of  target 
radius. 

The  stability  results  indicate  that  the  aspect  ratio  will  be 
limited  to  about  20  by  the  Rayleigh-Taylor  instability.  This 
implies  an  ablation  pressure  requirement  of  about  20  Mbars.  This 
result  can  be  obtained  from  a  simple  integration  of  the  equations 
of  motion  resulting  in  the  relation  "  1.5  <Pq>  (ARq/Rq). 

Thus  another  barrier  is  the  requirement  that  >  26  Mbars.  The 
region  of  acceptable  symmetry  and  stability  lies  to  the  right  and 
above  these  two  boundaries. 

The  requirement  of  lew  fuel  preheat  from  plasma  instabilities 
places  a  limit  that  scales  like  IX^.  This  curve  will  put  a  limit 
to  the  intensity  and  wavelength  from  above.  Unfortunately  we  do 
not  yet  possess  enough  knowledge  to  predict  the  absolute  location 
of  this  a  curve  on  the  graph  for  large  scale  reactor  plasmas  and 
this  remains  an  uncertain  parameter. 

If  we  now  add  results  for  efficiency,  we  can  see  that  the 
highest  hydrodynamic  efficiency  is  obtained  at  the  intersection  of 
the  smoothing  and  ablation  pressure  boundaries.  If  we  go  to 
shorter  wavelengths  along  the  symmetry  boundary  the  plasma  insta¬ 
bility  parameter  IX^  increases,  and  the  hydrodynamic  efficiency 
goes  down.  Increasing  the  wavelength  along  the  fixed  pressure 
boundary  curve  also  increased  IX^  and  decreases  hydroc^ynamic 
efficiency.  The  window  for  operation,  if  it  exists,  should  lie 
between  1/2  micron  and  2,7  microns  and  above  10^“*  W/cm^  irradiance. 
We  note  that  if  a  pessimistic  value  of  the  flux  limit,  f=0.03,  is 
used  the  smoothing  requirement  drives  the  optimum  region  to  2-3 
microns  where  the  efficiency  will  be  even  lower. 

Several  suggestions  have  been  made  to  overcome  the  smallness 
of  the  I-X  window  of  operation.  We  have  suggested  broadband  lasers 
to  reduce  plasma  instabilities,^  Nuckolls  has  suggested  a  series  of 
tenporal  changes  in  laser  frequency  starting  at  longer  wavelengths 
for  smoothing  and  ending  with  shorter  wavelength  for  maximum 
coupling  efficiency.^®  McCroy  has  suggested  that  the  direct  laser 
uniformity  may  be  controlled  to  better  than  ±10<>.  Skupsky  has 
suggested  that  there  is  enhanced  smoothing  in  the  flux-limited 
regime. ^  ®  One  encouraging  fact  is  that  absorption  by  inverse 


bremsstralung  is  at  lower  densities  in  longer  scale  length  plasmas. 
If  the  absorption  is  at  sufficientiy  low  densities,  then  thermal 
smoothing  beyond  the  critical  surface  may  take  place.  For 
sufficienl^'  large  radius  pellets  we  find  that  the  increased  inverse 
bremsstrahlung  of  shorter  wavelength  lasers  indeed  significant!;' 
increases  the  absorption  radius.  In  Fig.  11  we  show  the  location 
of  the  critical  surface  and  the  50*  absorption  surfaces  as  function 
of  time  for  red  and  green  laser  light  for  a  5  mm  target  under 
similar  icplosion  conditions  (shown  at  the  bottom  of  the  figure.) 
For  green  light  half  the  energy  is  absorbed  at  a  distance  3  times 
the  separation  distance  of  the  critical  to  ablation  surfaces.  This 
enhanced  smoothing  may  relax  somewhat  the  constraints  on  thermal 
smoothing.  However  we  note  that  the  long  scale  lengths  required 
for  this  phenomena  may  also  decrease  the  absorption  due  to  enhanced 
S3S.  Insufficient  information  is  yet  available  for  these  long 
scale  length  plasmas  to  draw  any  definitive  conclusions. 


TIME  (nsac) 

Figure  11.  Fatio  of  separation  distance  of  absorption  region  from 
critical  surface  to  target  radii  as  a  function  of  time 
for  red  and  green  light.  Lower  figure  is  radius  of 
ablation  surface  as  a  function  of  time.  The  conditions 
were  picked  to  match  the  implosion  velocities. 


Ve  are  still  left  with  a  n'-mber  of  uncertainties.  First  our 
results  are  based  on  hydrodynamc  sLuulations ,  which  neglects 
possible  phenomena  such  as  fast  electron  transport,  filamentation, 
and  other  plasma  instabilities.  As  yet  no  experimental  verifica¬ 
tion  exists  for  the  reduced  Rayleigh-Tkylor  results  and  little  is 
experimentially  known  about  plasma  coupling  at  reactor  scale  length 
plasmas.  What  our  results  do  indicate  is  that  unless  an  alterna¬ 
tive  is  found  to  smoothing  of  laser  nonuniformities,  short  wave¬ 
length  lasers  will  have  difficulties  providing  sufficient  smoothing 
for  reactor  sice  ICF  targets. 
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